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We demonstrate that the formalism for the calculation of the Jacobi last multiplier 
for a one-degree-of-freedom system extends naturally to systems of more than one degree 
of freedom thereby extending results of Whittaker dating from more than a century ago 
and Rao (Proc Ben Math Soc 2 (1940) 53-59) dating from almost seventy years ago. We 
illustrate the theory with an application taken from the theory of coupled oscillators. We 
indicate how many Lagrangians can be obtained for such a system. 
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1 Introduction 

Jacobi's Last Multiplier is a solution of the linear partial differential equation [3], HI [5j [6] , 
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where Y^iLi a i9 Xi is the vector field of the set of first-order ordinary differential equations for the 
iV dependent variables x i} namely 

A f = ^2a i (x 1 , . . . ,x n ) — = (2) 

or its equivalent associated Lagrange's system 

d^i dx 2 dx n 



(3) 



The ratio of any two multipliers is a first integral of the system of first-order differential equations 
and in the case that this system is derived from the Lagrangian of a one-degree-of-freedom system 
one has that [SI E] 

ffiT 

W = h ' L (4) 

Consequently a knowledge of the multipliers of a system enables one to construct a number of 
Lagrangians of that system. 



We recall that Lie's method [3 [13] for the calculation of the Jacobi Last Multiplier is firstly 
to find the value of 

r e» 



A = det 



(5) 



in which the matrix is square with the elements being the vector field of the set of first-order 
differential equations by which the system is described and the elements, Sy, being the coefficient 
functions of the number of symmetries of the given system necessary to make the matrix square. 
If A is not zero, the corresponding multiplier is M = A -1 . There is an obvious corollary to 
the results of Jacobi mentioned above. In the case that there exists a constant multiplier, the 
determinant is a first integral. This result is potentially very useful in the search for first integrals 
of systems of ordinary differential equations. 



The relationship between the Jacobi Last Multiplier and the Lagrangian for a one-degree- 
of-freedom system is perhaps not widely known although it is certainly not unknown as can 
be seen from the bibliography in [16J. Given a knowledge of a multiplier, (jSJ) gives a simple 
recipe for the generation of a Lagrangian. The only possible difficulty is the performance of the 
double quadrature. In this paper we consider that the vector fields of the system of equations 
and symmetries are known and that we seek the multiplier. From another direction one could 
know the multiplier and all but one of the symmetries. From (jSJ) the remaining symmetry can be 
determined [15J. A knowledge of the multipliers of a system enables one to construct a number 
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of Lagrangians of that system. Considering the dual nature of the Jacobi Last Multiplier as 
providing a means to determine both Lagrangians and integrals one is surprised that it has not 
attracted more attention over the more than one and a half centuries since its introduction. The 
bibliography of [16J gives a fair indication of its significant applications in the past to which we 
may add HU, HD|, [TBI, jig, TO PU- 



2 Theoretical Development 

In Rao [22] there is a development of the work of Whittaker [6] [pp 276-286] on the theory of the 
Last Multiplier of Jacobi for multidimensional systems. The discussion by Whittaker is limited to 
dynamical systems without a potential. The work of Rao is restricted to dynamical systems with 
a Lagrangian which is a quadratic form in the generalised velocities. It does provide a clearer 
path to an understanding of the methodology than that provided by Whittaker. In particular 
Rao gives specific formulae for the equations to be satisfied by the various multipliers. We see 
in the sequel that a considerable simplification occurs, but we pay due acknowledgement to the 
inspiration afforded by that paper. Although the theory which we develop here is applicable to 
dynamical systems of all dimensions, we confine our attention to two dimensions to keep the 
mathematics as simple as possible so that the essential ideas are quite evident. We assume that 
we have a Lagrangian, L(t, qi,q2, Qi, fa), in a standard notation. We write the simplified form of 
the corresponding Euler-Lagrange equations as 

Qi = fi(t,qi,q 2 ), (6) 
fa = f2(t,qi,q 2 ), (7) 

where we note that the generalised force is independent of the generalised velocities. The corre- 
sponding vector field is 

X l = d t + fad qi + q 2 d q2 + fxd qi + f 2 d ia . (8) 
The actual Euler-Lagrange equations are 

d 2 L d 2 L 3 2 L 3 2 L d 2 L „ dL 

+ + 77^^2 + -TTTjfl + ^T^/ 2 Tj = (9) 



dq\dt dqidqi dq\dq2 dqj dq\dq2 dqi 
8 2 L d 2 L d 2 L d 2 L „ d 2 L „ 8L 

dq 2 ot oq 2 oq 1 oq 2 oq2 oqidq 2 aq{ Oq 2 

in which we have replaced q\ and fa by the right sides of the equations of motion above. 
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We follow Rao [22] in defining the connection between the last multiplier and the Lagrangian 

as 

ie, we assume the same relationship as for the one-degree-of- freedom case so that, once the 
multiplier has been calculated, the Lagrangian follows by a double quadrature. 

We differentiate (jUJ) and ffTUl) once with respect to both q\ and q 2 . We illustrate the calculation 
in the case of the differentiation of (Q with respect to q 2 . We obtain 

d 3 L d 3 L . 8 3 L 

dqi dq 2 dt dq x dq 2 dqi 1 dq\ dq 2 dq 2 ^ 

| d 2 L + 8 3 L + 8 3 L j _ d 2 L _ Q 
dq\dq 2 dq 2 dq 2 dq\dq 2 dq 2 dqi 

in which we have used the independence of f\ and f 2 from q 2 (in this instance). We use the 
definition, (JTT1) . for M\ 2 . Then (fT2|l becomes 



| (Ml2 ) + A (M 12(il ) + A (Ml2 j 2 ) + A (M 12/l ) 

In the cases of M u and M 2 2 the ultimate and penultimate terms in the equations corresponding 
to f|T3|) cancel. The other equation for M 12 has the subscripts reversed so that, when the two 
equations are added, the terms vanish. Consequently each multiplier is a solution of the equation 

m < M > + £ < M «>> + w 2 < M *> + 4 {Mh) + 4 (M/2) = °' (14) 

Allowing for the variations in notation (JT3J) is the same equation as (HI), ze, the number of degrees 
of freedom of the system only affects the number of terms in the equation. 



3 Example: A two-dimensional coupled oscillator 

3.1 Equivalent Lagrangians 

The collinear motion of two particles of equal mass connected to each other by a spring and 
to fixed points by two other springs of equal spring constants is modelled by the fourth-order 
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system^. 



qx = - Ql{qi - g 2 ) (15) 

g 2 = - g 2 ) - ^??2- (16) 

It is evident that a solution of (fl~4l is simply M = constant. We take the constant to be unity 
for Mn and M 22 and zero for M i2 . Then the corresponding Lagrangian is 

L=\(q\ + qf) + hi(t, q 1: q 2 )qi + h 2 (t, g x , q 2 )q 2 + h 3 (t, qi,q 2 ), (17) 

where the three functions, hi, % — 1, 3, are functions of integration. When we substitute L, (fTT|) . 
into the general Euler-Lagrange equations, we require that we obtain equations ffTB^) and ffTB]) . 
This requirement imposes the constraints 



dh dh 2 Q 



<9g 2 dqi 

- ti( qi - n\{ qi - q 2 ) = (19) 



^1 ^3 rV 2 ,-,2 



<9t <9gi 

+ ^(g 1 -g 2 )-n 2 g 2 = 0, (20) 



^2 9^3 n2l . , n2 



<9t <9g 2 

where the constraint (ITS]) is common to both equations. The constraint (ITS]) is the two-dimensional 
version of the condition for the functions to be the divergence of some arbitrary function, 
g(t,qi,q 2 ), ie 

hi = and h 2 = (21) 

dqi oq 2 

When ( I2TI) is substituted into (jT9j) and (I20l . one obtains 

<>>r/ <J/ ' 3 -o 2 g 1 -Q 2 : (g 1 -g 2 ) = o (22) 



<9t<9g 2 <9g 2 
which can be consistently integrated to give 



dtdqi dqi 

+ ^(gi-g 2 )-Q^ 2 = 0, (23) 



^ _ 1 

at 2 



+ ^2) (9? + ?2 2 ) - 20^ig 2 ] + s(t), (24) 



1 A simplified version of this system can be found in [Hj[pp 479-481] 
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where s(t) is an arbitrary function which may as well be set at zero, since it makes no contribution 
to the Euler-Lagrange equation, or included in the other arbitrary function g(t, qi, q 2 ) which also 
makes no contribution to the Euler-Lagrange equation. 



Eventually we have arrived at the Lagrangian 



L = | (<g + it) + 9 - \ [(«? + (<Z? + ql) - 2fifo<& 



(25) 



ie the arbitrariness in the Lagrangian can be expressed as a total time derivative. Such a 
Lagrangian has been termed 'gauge variant' [23] and is notable in that the presence of the 
arbitrary function has no effect upon the number of Noether point symmetries [20] . In this 
respect it could be regarded as part of the boundary term in the way Noether put it in her 
formulation of her theorem [23]. The class of Lagrangians described by (1231) is an equivalence 
class. 



3.2 A Plethora of Multipliers 

Jacobi proposed to obtain the last multiplier from the solution of the first-order linear partial 
differential equation, (fUJ) . In 1874 Lie [7] showed that one could use point symmetries to 
determine last multipliers. Lie's method for the calculation of the Jacobi Last Multiplier is 
firstly to find the value of the determinant in (|14p . The number of symmetries required depends 
upon the number of elements in the vector field, which supplies the first row of the matrix, since 
the matrix must be square. Here we consider that the vector fields of the system of equations 
and symmetries are known and that we seek the multiplier. The Lie point symmetries of system 
(USHISD are 

Ti = q 2 d qi + qid g2 + q 2 d qi + qidg 2 

T 2 = exp [Qit] [(d qi - d q2 ) + iQ (d qi - d q2 )) 

r 3 = d t 

T 4 = qid qi + q 2 d g2 + qid qi + q 2 d q2 

T 5 = exp [-Qit] (d qi - d q2 - iQ (d qi - d q2 )] 

T 6 = exp [iQxt] [(d qi + d q2 ) + iQ x (d qi + d q2 )} 

T 7 = exp [-iSl x t] [(d qi + d q2 ) - iVt^t (d qi + d q2 )} , 

2 At that time the only symmetries considered were point symmetries. With the extension of the concept of 
infinitesimal transformations generated by symmetries depending upon derivatives and integrals there has been 
a natural extension of Lie's valuable contribution to the subject of the Jacobi Last Multiplier. 
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where Q = y Clf + 2tt 2 , which is the minimal number of Lie point symmetries for an autonomous 
second-order two-dimensional linear system [25] - Despite the fact that the system possesses 
only the minimal number of Lie point symmetries there are thirty-five different determinants 
to be evaluated. Naturally not all of these can be expected to be nonzero or different, but one 
can easily imagine that this would be the case for a substantial proportion since in the case of a 
single linear oscillator the number of different multipliers is fourteen out of a possible twenty-eight 
determinants 



It is not our intention to display a large number of multipliers. We simply wish to illustrate 
two points. The first is the determination of a Lagrangian for the system ffl5| /fT6|) differing from 
the class represented by (125]) which is of a standard form. The second is to display a first integral 
for the system by taking the ratio of two autonomous multipliers. From the Lagrangian (j2"5l) 
one can construct an autonomous Hamiltonian which is a first integral. Thus we have two first 
integrals and, as if we did not already know it, the system is integrable from the Theorem of 
Liouville. 



If we take the symmetries Ti, r 2 , T 3 and T 5 , the matrix for which the determinant is to be 
evaluated is 



C 



1235 



1 q x q 2 - ^1(<?1 - 92) £M<Zl - 92) - Qlq 2 

g 2 qi <?2 qi 

exp [Qit] — exp [Qit] iQ exp [Qit] —iQ exp [Qit] 

1 

exp [— Qit] —exp [—Qit] —iQ exp [—Qit] iQ exp [—Qit] 



We obtain the multiplier 



(26) 



JLM 1235 = — 



2iQ [(qi + q 2 ) + fi? (?i + <&) 
Similarly by taking the symmetries Ti, Fq, r 3 and Tj we obtain 

JLM lw ' 



(27) 



Obviously an integral is 



2VL 2 {q x - q 2 f + Vt x {q x - q 2 f + {qx - <j 2 ) 2 



2iQ {qx + q 2 f + (<?i + q 2 ) 



(28) 



1637/1235 



2tt 2 {qx - q 2 f + fix {qx - q 2 f + {qx - q 2 f 



(29) 
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We can obtain a Lagrangian from JLM1235, say, by integrating twice with respect to q\. We 
have 



1 



J 1235 



2iQVtx (q 1 + q 2 



(91+92) arctan 



4zQ 



flj (9i + 9 2 ) 2 + (9i + 92) 

n? (91 + 9 2 ) 2 



9i + 92 \ 
Pi (91 + 92)/ 

f Rt(t, qt,q 2 , q 2 )qi + Rzit, 9i, 92, 92) 



(30) 



where R\ and R 2 are functions of integration. When we apply the Euler-Lagrangian condition 
to ( 130]) and substitute the required equations of motion, (!T5]/ fl6l). we obtain the conditions 



2Q(9i + 9 2 ; 



+ 



^2 (9i - 92) - ^l92 



dR x . dRi dR x dR 2 n 
+ q%— — = 



and 



2Q(9i + 92) 
d 2 R 1 



Sl\qx + ^2 (9i - 92) 



dq 2 dq 2 dt dq 1 
dR x . dR x dR 2 



(31) 



+9i92 



dq 2 dq<2 



+ 9r 



2 d 2 R! . d 2 R ± 



dqidqi 



+ 9l- 



dt 2 



dq 2 
= 0. 



9i 



dqi dq 2 



+ q x Vt 2 (qi - q 2 ) - n x q 2 



d 2 R x 
dql 



(32) 



This second equation is a quadratic in qi which is not one of the variables in Ri and R 2 . The 
three coefficients give 



dqi 



+ 



^2 (9i - 92) - ^l92 



2<5(9i + 92) 



d 2 R 1 
dqtdq 2 

d 2 Ri d 2 R x d 2 Ri 
dq 2 dq 2 dq 2 dt 2 
dR x dR 2 



Q^i + Vt 2 (gi - q 2 



dq 2 dq 2 





0. 



(33) 
(34) 
(35) 



We note that R 2 appears only in (I3TI) and (I35I) and in a very symmetrical way at that. Equations 
(|33|) and (I34|) are easily satisfied by setting Ri equal to zero. It then follows that we can take 



R<2 = 2Q l0g ^ + ^ ' 



(36) 
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4 Discussion 



The last multiplier of Jacobi has in the main been sadly neglected. It is evident from the slight 
literature [HI [22] that its application to systems of more than one degree of freedom was either 
ignored or misunderstood. In §2 we saw that there was no difficulty to derive the basic equation 
for the multiplier. It is the same equation as for the determination of a multiplier for a one- 
degree-of-freedom system. In the context of our presentation the main implication of that result 
was that Lie's method was equally applicable to systems of more than one degree of freedom. For 
the sake of simplicity of presentation we confined our attention to just two dependent variables, 
but it is quite obvious that the result is general. When we considered the specific example of 
the system, (fT5]4TB|) . we used obvious solutions of (jHJ) to obtain a class of Lagrangians which 
includes the usual Lagrangian for such a system. Furthermore we demonstrated the use of the 
symmetries of the system to construct two other multipliers and hence a second integral of the 
system^. Consequently the system, (fl5|/fT6l). is integrable in the sense of Liouville. This is to 
be expected given that the system is a pair of coupled linear oscillators. In addition to the 
'standard' Lagrangian^, (j25p . we explicitly demonstrated the existence of another Lagrangian. 
One infers the existence of many. 

It is amazing to think that in the case of a two-dimensional linear oscillator, which admits 
a 15-dimensional Lie point symmetry algebra (the maximum for a two-dimensional system), the 
possible Lagrangians have to be sought among 1365 different determinants! 

Further work is in progress in order to address the case of dissipative systems. 
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